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Abstract: As part of global climate change an accelerated hydrologic cycle (including an increase in 
heavy precipitation) is anticipated ([21], [22]). So, it is of great importance to be able to quantify high- 
impact hydrologic relationships, for example, the impact that an extreme precipitation (or temperature) 
in a location has on a surrounding region. Building on the Multivariate Extreme Value Theory we 
propose a contagion index and a stability index. The contagion index makes it possible to quantify the 
effect that an exceedance above a high threshold can have on a region. The stability index reflects the 
expected number of crossings of a high threshold in a region associated to a specific location i, given 
the occurrence of at least one crossing at that location. We will find some relations with well-known 
extremal dependence measures found in the literature, which will provide immediate estimators. For 
these estimators an application to the annual maxima precipitation in Portuguese regions is presented. 

1 Introduction 

The need to model and predict environmental extreme events such as hurricanes, floods, droughts, heat 
waves and other high impact events, which can lead to a devasting impacts, ranging from disturbances 
in ecosystems to economic impacts on society as well as loss of life, motivated the modeling of spatial 
extremes. 

A common method of modeling spatial extremes is through max-stable processes, which are an infinite 
dimensional generalization of multivariate extreme value distributions (Haan (1984) [4], Vatan (1985) [23] . 
de Haan and Pickands (1986) [3]). 

Max-stable processes can be, for example, good approximations for annual maxima of daily spatial 
rainfall (Smith [18] , Coles [2], Schlather |14| . among others) and therefore have been widely applied to 
real data. 

Briefly, a max-stable process X = {A"i} igKd is the limit process of maxima of i.i.d. random fields 
i j ^ 1- Namely, for suitable \a n 

(i)> 0)n>i an d {^n(i)} n >i sequences of real constants, 



provided the limit exists. 

We shall consider d = 2, that is X = {^i}i G iR2 • The distribution of . . . ,X- lk ) is a Multivariate 
Extreme Value (MEV) distribution Ga, A = {ii, . . . , i*,}, and we can assume, without loss of generality, 
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that the margins of X have a unit Frechet distribution, F(x) = exp(— x 1 ), x > (Resnick [12]). The 
distribution Ga can then be defined by 



where Va denotes the exponent function of the MEV distribution Ga- 

The exponent function summarizes the extremal dependence structure of Ga and the scalar Va(1, •••,!) 
defines the extremal coefficient ea detailed in Schlather and Tawn (2003. [To]), which summarizes the ex- 
tremal dependence between the variables indexed in the region A. This coefficient takes values between 

1 and k, with a value of 1 corresponding to complete dependence and a value of k corresponding to 
complete independence. Its value can be thought of as the number of effectively independent locations 
among the k under consideration. 

If we consider A = {i, j} we find the extremal coefficient of Tiago de Oliveira (1962/63, [20]) which 
is related with the bivariate upper tail dependence coefficient A{;j} = lim u -|-i P (F(Xi) > u\F(Xj) > u), 
introduced in Sibuya (1960, [IT]), as e{;j} = 2 — A{; n. 

Although these measures are very useful to analyze the dependence among extremal events, there 
remain important questions to be answered, for example, the influence of an extreme event on the 
regional smoothness of X and the contagion effect of an extreme event at a specific location over the 
variables of X indexed in a region of M. 2 . 

In Section 2 we propose a measure, called contagion index, that enables to quantify the impact that 
an exceedance of a high threshold can have on a region and we present its relation with bivariate extremal 
coefficients. 

Clearly an extreme event could affect the smoothness of a random field over a region so, in Section 
3 we propose a stability index on a region A associated to a specific location i, defined as the expected 
number of crossings of a high threshold u in A associated with i, given that there is at least one crossing 
in A associated to i. We also present some properties of this coefficient. 

Section 4 is devoted to illustrate the previous measures in a max-stable M4 random field. 

Based on relations of our indices with well-known dependence measures, for which estimators and 
respective properties have already been studied in the literature, in Section 5 we present estimators for 
the stability and contagion indices. We end with an application to the annual maxima precipitation in 
Portuguese regions. 

2 Contagion Index 

The occurrence of an extreme event at a given location i may spread throughout a region of locations. 
In this section we define a measure for assessing the effect of an exceedance above a high threshold u at 
a specific location i on a region A of locations. 

Definition 2.1 Let X = {Xi} ieR 2 be a max-stable random field with unit Frechet margins F and A a 
region of M 2 . The contagion index from the location i to the region A is defined as 





(1) 
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The CI(A, i) is the conditional expected number of exceedances above a high threshold u in A, given 
X\ exceeds u, that is, the CI(A, i) measures the impact that the event {X{ > u} has on the region A. 

We remark that the conditioning location i does not necessarily have to be in the region A. 

The following proposition states that CI (A., i) relates with the bivariate extremal dependence 
coefficients j G A, and the tail dependence coefficients Ajjj}, j G A. 

Proposition 2.1 For any max-stable random field with unit Frechet margins F, i G M 2 and A C R 2 , we 

have 

CI(A,i) = E A {y} = 2|A|-5> {y} . 
jeA jeA 

Proof: Observe that 

C/(A,i) = ^limP(F(X j )>u|F(X i )>u) 

= E x m = E ( 2 " = 2 1 A l - E e iiJ}- 

jeA jeA jeA 

□ 

An C/(A,i) close to |A| means that i has a high influence on A, while an CJ(A, i) close to zero 
implies a negligible influence of i on A. In other words, the higher the index, the higher the contagion 
effect of the event {X{ > u} on the region A. 

Remark 1 To gain some intuition for this measure, as a device for measuring dependence, consider two 
polar cases: 

• Case 1. If X\ is independent of Xj, for each jGA, then CI(A, i) = 0. 

• Case 2. If, for each j G A, X- } and X\ are totally dependent, then CJ(A,i) = |A|. 



Remark 2 We can extend the CI(A, i) to the contagion index from a region A to a region B, as follows 



C/(A,B)=lim£PTlI m x jW 
VjeA 



]J{F(X i )>u} 



ieB 



This measure is related with the multivariate upper tail dependence coefficient (Schmidt (2002, |13|): Li 
(2009, [IE]) ; Ferreira (2011, [6])), defined as 



A A , B = hmP f|{F(Xj)> U } 



nw)>«} ) , 

ieB 



in the following way 



CI(A,B) = J2 
jeA 



E ^ J c B (-l) |J|+1 A J , { j } 

£B 



When we take A = B, we obtain the fragility index (FI) of the region A. The FI was introduced in 
Geluk et al. (2007, [9]) to measure the stability of a stochastic system. The system is called stable if 
FI = 1, otherwise it is called fragile. 
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3 Stability index of a region A 



In order to analyze the regional smoothness of a random field associated to a specific location we 
propose the following measure. 

Definition 3.1 Let X = {X^}-^^ be a max-stable random field with unit Frechet margins F and A a 
region ofM 2 . The stability index of the region A associated to a specific location i £ M?, SL(A,i), is 
defined as 



5/(A,i) = \imE \J2 I{FW)<«<^(Aj)} 
W 



2 ^{F(Xi)<u<F(Xj)} > 
J6A 



The SI( A, i) is the conditional expected number of crossings (above a high threshold u) in A from a 
specific location i, given that there is at least one crossing in A from i. 

If a max-stable random field X does not vary smoothly over a region A, we will expect a large number 
of crossings of a high threshold in A associated to a specific location i. A higher number of crossings 
signifies increased instability. 

The next results highlight the connections between SI( A, i) and the extremal coefficients. 

Proposition 3.1 For any max-stable random field with unit Frechet margins, i € IR 2 and A C M?, we 

have 

s/(A , i) = 5^>-l A l-l A l-«(A,i) 



-{i}UA e {i}UA 



Proof: Since 

E Yl 2 {^(*i)<«<F(Xj)} 

VjeA 



J6A 
it follows that 

I A| u - £ jeA « e < IJ} _ E jeA - I A 



5/(A,i) = lim 



-{i}UA 



□ 



Proposition 3.2 Under the conditions of Proposition 3.1., we have 

Y^jeAtm ~ |A| < SI{A/i) < EjeAtm - IAI 



|A| + 1 " ' ' Vj. A' ;i..i; 

Proof: Just observe that 6{i}uA < |A| + 1 and e|i} UA > VjeA e {ij}- ^ 
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4 An M4 random field 

In this section we derive the contagion and stability indices for a particular class of max-stable processes 
known as multivariate maxima of moving maxima, or M4 processes for short (Smith and Weissman 
(1996), [IH])) which are particularly well applicable in a time series or spatial process context. An M4 
random field X = {Xi} ieZ2 is defined as 

Xi = max max ai mi Zi i_ m , i G Z 2 , (3) 

1>1 -oo<m<+oo 



where {Zz n}z>i,nez is a family of independent unit Frechet random variables and, for each i € Z , 

+oo +oo 

{aimi}i>i,mez are non-negative constants such that a\ m \ = 1. By considering that the distribu- 

1=1 m=—oo 

tion of (Xij , . . . , X\ k ) is characterized by the copula 

+oo +oo 

C(u h ,...,u ik ) = H 11 /\ u?™, Ui, € [0,1], j = l,...,k, (4) 

1=1 m=-co ie{ii ,...,ifc} 

it was shown in Fonseca et al. (2012, [S]) that the random field X = {Ai} igz2 is max-stable and the 
exponent function of the distribution of (X^ , . . . , X\ k ) is given by 



+oo +oo k 

VA(xi,...,x k ) = ^2 \/ (^j 1 ^!^! , Xj eM, j = l,...,k, A = {i!,...,i k }. 

1=1 m=—ooj=l 

So 



-oo +oo 



CI(A,i) = 2|A| - Y^Yl Yl ( a irniVa lmi ) 

jeA 1=1 m=-oo 

all< ^ + + 

SI(A i) — ^j eA ^ =1 ^m=-oo( a lrni V a; m j) — |A| 



E/=T Em=-oo (VjeA a lmj V a h 

To illustrate the computation of the contagion and stability indices we shall consider, in what 
follows, examples with a finite number of signature patterns (1 < I < L) and a finite range of sequential 
dependencies (Mi < m < M 2 ). 

Example 4.1 Lets consider that for each location i € Z 2 with even abscissa we have 
«iii = f> oi2i = \ and otherwise am = | = 1 — ai 2 \- The values of (ani,ai2i) determine the 
moving pattern or signature pattern of the random field, which in this case corresponds to one pat- 
tern (L = 1). 

Let A = {sj(i) : i = (3,3) A j G {1,2,..., 8}}, where Sj(i), j = 1, . . . , 8, denote the neighbors of i 
defined as follows: 

si(i) = (ii + l,i 2 ) s 2 (i) = i+l s 3 (i) = (h, i 2 + 1) s 4 (i) = (h - 1, i 2 + 1) 
s 5 (i) = (it - l,i 2 ) s 6 (i) = i-l s 7 (i) = i 2 - 1) s 8 (i) = (ii + l,i 2 - 1). 
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The matrix of the bivariate extremal coefficients, e 



, e {sj(i),i}j J 



, 8, provide insight into the likelihood 



of contagion from i = (3, 3) to its neighbors Sj(i), j = 1, . . . , 8, although without specifying the size of 
the impact which is given by CI (A, i). We obtain 



£ {s4(i),i} 


e{s 3 (i),i} Hs 2 (i),i} 




31 
20 


1 


31 
20 


£ {s5(i),i} 






31 
20 


1 


31 
20 


e {«6(i),i} 


e {sr(i),i} e {^8(i),i} . 




31 
20 


1 


31 
20 



and 

C7(A,i)=4,7. 

The stability index of the region A associated with i = (3, 3) is given by 



Example 4.2 Now, we shall consider one example with two signature patterns (L = 2 

Lets assume that for each location i € 1? we have am 
i 



Ol2i 



a 13i — 5) 

«2ii = «22i = j^j , a23i = k if both coordinates are odd and am = \,a\2i = «i3i = g, 
a 2ii = «22i = «23i = | otherwise. The values of (am, ai2h a i3i) an d («2ii 5 fl22i; «23i) define the two 
signature patterns of the random field. 

For A = {(2, 4), (3, 4), (4, 4), (5, 4)} and i = (3, 3) we obtain C7(A, i) = f§ and SI(A, i) = ff . 



We next focus on the estimation of the stability and contagion indices. 



5 Estimation 

As previously stated, the contagion and stability indices relate with the extremal coefficients of Tiago de 
Oliveira (1962/62, [20]) and Schlather and Tawn (2003, [E]), which can be expressed by the exponent 
function given in (1). 

There are several references in literature on the estimation of the exponent function. For a survey we 
refer to Krajina (2010, [10]) and Beirlant (2004, [I]). 
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Figure 2: Simulation of the M4 random field with L = 2. 1 < m < 2 and for constants defined in Example 4.2. 



It is known that parametric estimation methods are efficient if the distribution model under consid- 
eration is true, but they suffer from biased estimates otherwise. Non parametric estimation procedures 
avoid this type of model error. However, they are usually based on an arbitrarily chosen parameter k 
corresponding to the number of top order statistics to be used on the estimation of a high quantile of 
F, which relates to the usual variance-bias problem: if k is too small, then the estimator tends to have 
a large variance, where if k is too large, then the bias tends to dominate. Some methods of choosing an 
optimal k are discussed in Einmahl et al. (2006, [3]). 

In order to overcome the problem of the optimal choice of k, Ferreira and Ferreira (2011, |7j) developed 
another approach. Based on the following relation 

eA = y A (l,l,...,l) = I j^^y , where M(A) = \J F^X,), 
the estimator of ca proposed in Ferreira and Ferreira (2011, [7]) is defined as 



M(A) 



1 - M(A) 



where M(A) is the sample mean, 

i n 

^) = TV^ (J) : 



n 

3=1 ieA 



and Fi, i € A, is the (modified) empirical distribution function of F{, 



1 n 

3=1 



With this estimator of the extremal coefficient and the relations established in Propositions 2.1 and 
3.1 we propose, respectively, the following estimators for the contagion index C/(A, i) and the stability 
index SI (A, i), 

C/(A,i)=2|A|-^? {iJ} 
jeA 

and 
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e {i}UA 

which are consistent given the consistency of the estimators eim and e^ju^ already stated in Ferreira 
and Ferreira (2011, [7]). 

The results of the application of these estimators to the Examples 4.1 and 4.2 are presented in the 
following tables: 





CI 


CI 


MSE 


Example 4.1 


4.7 


4.70806 


0.01259 


Example 4.2 


3.2667 


3.26684 


0.0006 



Table 1: Results with 100 replications of 1000 i.i.d. max-stable M4 random fields of the Examples 4.1 
and 4.2 where CI denotes the true values of the contagion index, CI the estimated values and MSE the 
estimated mean squared error. 





SI 


SI 


MSE 


Example 4.1 


2.12903 


2.1222 


0.00234 


Example 4.2 


0.61972 


0.62029 


0.00034 



Table 2: Results with 100 replications of 1000 i.i.d. max-stable M4 random fields of the Examples 4.1 
and 4.2 where SI denotes the true values of the stability index, SI the estimated values and MSE the 
estimated mean squared error. 

As we can see the estimated values are very close to the true values of the coefficients. These results 
show that the simple non-parametric estimators CI and SI are a promising tool for assessing regional 
contagion effects and regional smoothness for these random fields. 

6 An application to precipitation data 

We now illustrate the estimation of the contagion and stability indices through an application to 
annual maxima values of daily maxima precipitation recorded over 32 years, in six Portuguese stations 
(Figure 3), provided by the Portuguese National System of Water Resources (http:/ /snir h.pt[ ). 

Since the data are maxima over a long period of time, we assumed that they are independent over 
the years in each location. We also assumed that the random field is max-stable with unknown marginal 
distributions so data were previously transformed at each site so that they have a standard Frechet 
distribution. 

The estimated values of the contagion and stability indices from Lagoa Comprida (located at the 
highest altitude) to the regions {Gouveia, Oliveira do Hospital, Seia} and {Penamacor, Barragem Cabego 
Monteiro} are presented in Table 3. 

The results suggest that Lagoa Comprida has a higher influence on the region {Gouveia, Oliveira do 



8 




200000 220000 240000 260000 280000 300000 



Figure 3: The locations of the stations where precipitation data were collected, obtained from Portuguese National 
System of Water Resources (left) and their representation in Lambert coordinates (right). 



A 


CI (A, Lagoa Comprida) 


SI (A, Lagoa Comprida) 


{Gouveia, Oliveira do Hospital, Seia} 


0.96688 


0.8959 


{Penamacor, Barragem Cabego Monteiro} 


0.00887 


0.75741 



Table 3: Estimates of the contagion and stability indices from Lagoa Comprida to the regions {Gouveia, 
Oliveira do Hospital, Seia} and {Penamacor, Barragem Cabego Monteiro}. 



Hospital, Seia} in terms of precipitation amounts and this region is smoother when compared to region 
{Penamacor, Barragem Cabego Monteiro}. 
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